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Inclusive charm meson production cross sections in the deep inelastic scattering and photoproduction regimes 
are compared with QCD leading and next-to-leading order (NLO) calculations. The NLO predictions are signifi- 
cantly below the data in some parts of the measured kinematic range. Angular distributions of dijet events with 
charm show clear evidence for the existence of charm originating from the incoming photon. The charm frag- 
mentation function is measured for the first time at HERA. Various fragmentation ratios and the fragmentation 
fractions of the low-lying charm hadrons are determined and compared to previous e + e~ results. 



1. Introduction 

During 1996-2000 HERA collided electrons or 
positrons (E e =27.6 GeV) with protons (E p be- 
tween 820 and 920 GeV). Open charm (c) pro- 
duction, which has been extensively studied, oc- 
curs in few steps: a) Hard process, e.g. boson- 
gluon fusion (BGF), where a photon (7) or Z 
boson emitted from the incoming electron fuses 
with a gluon (g) from the proton, producing 
a cc pair; b) Initial/final state parton shower 
development; c) Fragmentation of a final-state 
parton into a hadron. Two kinematic regions 
have been explored: 1) Deep inelastic scatter- 
ing (DIS) with photon virtuality Q 2 > 1 GeV 2 , 
where the scattered electron is visible in the 
main detector; 2) Photoproduction (PHP) with 
< Q 2 >« 3 • 10 -4 GeV 2 , where the virtual pho- 
ton radiated from the incoming electron is quasi- 
real. 

The large mass of the c quark provides a "hard" 
scale needed for the comparison of data to QCD 
predictions. In leading order (LO) QCD, two 
types of processes are responsible for charm PHP: 
Direct photon processes, where the photon inter- 
acts as a point-like particle and resolved photon 
processes, where the photon acts as a source of 
partons. BGF is the dominant direct photon pro- 
cess. Charm quarks present in the parton distri- 



butions of the photon lead to LO resolved charm 
excitation processes like eg — > eg. 

Various NLO calculations exist for comparison 
with charm HERA data: 1) fixed-order (FO) ap- 
proach PP, where only u,d,s are active flavours 
in the photon and proton; 2) resummed (RS) 
approach |2I,'{| . where a "massless" c quark is 
also an active flavour; 3) matched (FONLL) cal- 
culation 4 , which incorporates mass effects up 
to NLO and the resummation of pt logarithms 
up to next-to-leading logarithm (NLL) level. All 
these approaches are based on the DGLAP evolu- 
tion 5 j. The CASCADE Monte Carlo (MC), 
based on the CCFM evolution 0, is a more recent 
tool to compare with the data. 

Charm can be tagged by reconstruct- 
ing charm mesons via decays such as 
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D*+ -> D°n+ -> (!f-7r+)7r+, D+ -> K-%+w+, 
D° -> K-W+, D+ -> cV+ -> (K-K+)ir+, 
where charge conjugate states are included, or by 
measuring semileptonic electrons or muons from 
charm. Reconstructing charm hadrons via their 
decay vertices with a silicon tracker significantly 
reduces the background. 



2. Inclusive charm meson production 

Inclusive production of the charm mesons D* + , 
D°, D s and D+ in DIS has been studied |Hj 
with the HI detector, using a sample correspond- 
ing to an integrated luminosity of 48 pb . A 
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silicon tracker separates the production and de- 
cay vertices of the pseudoscalar mesons. Dif- 
ferential cross sections have been measured for 
the D + and D° mesons and compared with the 
AROMA LO MC predictions in the kinematic 
range 2 < Q 2 < 100 GeV 2 , 0.05 < y < 0.7, 
PtiP) > 2.5 GeV and \ij(D)\ < 1.5, where y is 
the electron inelasticity and pt and r] are, respec- 
tively, the transverse momentum and pseudora- 
pidity of the D meson. r\ is defined to be pos- 
itive in the proton beam direction. Results for 
the D + channel are shown in Fig.l. The lower 
shaded bands indicate the AROMA beauty con- 
tribution, scaled by the excess factor of data over 
MC, as found in the HI publication JU|. The 
distributions are well described by the LO QCD 
simulation, both in normalisation and shape. 
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Figure 1. Differential cross sections for D + as a 
function of the D + (left) and event (right) vari- 
ables. The AROMA MC predictions are given by 
the shaded bands. 



Inclusive PHP of D* + mesons has been mea- 
sured JT] with the ZEUS detector in the kine- 
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Figure 2. Differential cross sections for inclusive 
D* production as a function of px, f], W and 
z. FO predictions with nominal parameters are 
given by solid histograms. Upper (lower) dashed 
histograms correspond to upper (lower) bounds 
of the predictions. NLL predictions are shown 
by solid curves and shaded bands (dash-dotted 
lines) for the AFG (GRV) photon structure func- 
tion parametrisation. Direct photon NLL predic- 
tions are given by the dotted lines. 



1.9 < p t (D*) < 20 GeV and \rj(D*)\ < 1.6, us- 
ing an integrated luminosity of 79 pb" 1 . The mea- 
sured differential cross sections were compared 
with the NLO predictions FO (TJ, RS 2 and 
FONLL gj. In fig.2 the distributions da/dp T , 
da/drj, da/dW and do jdz are compared with the 
FO and RS calculations, where z(D*) is the frac- 
tion of the photon energy carried by the D* me- 
son in the proton rest frame. Theoretical uncer- 
tainties obtained by varying the charm mass and 
renormalisation scale are large, in particular for 
the NLL predictions. The central FO predictions 
are below the data, mainly for r\ > and low 
z. The NLL calculations are closer to the data. 
In particular NLL is better than FO for da jdz 
and for r\ > 0. A significant resolved contribution 
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is required in the NLL predictions, which show 
some sensitivity to the photon structure function 
parametrisation. 

The precise ZEUS data enable measurements of 
double differential cross sections. In Fig. 3 the rj 
distributions for four regions in px are compared 
with FO and FONLL predictions. The data is 
close to the upper band of the predictions and is 
significantly above the FO and FONLL calcula- 
tions at medium px and positive rj, where even 
the upper bounds are below the data. FONLL 
predictions are close to the FO ones for low px, 
but below FO for large px, where FONLL should 
do better. 
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Figure 3. Differential cross sections da / drj for in- 
clusive D* production for four px regions. FO 
predictions with nominal parameters are given 
by solid histograms. Upper (lower) dashed his- 
tograms correspond to upper (lower) bounds of 
the predictions. FONLL predictions are shown 
as dotted curves with uncertainties given by the 
shaded bands. 



3. Charm dijet angular distributions 

Charm-dijet PHP events enable the study of 
the photon structure, in particular it's charm 
content. Inclusive PHP of D*-dijet events has 
been measured with the ZEUS detector 

in the kinematic region Q 2 < 1 GeV 2 , 
130 < W < 280 GeV, E 3 T et > 5 GeV, 
\rf et \ < 2.4, M j:j > 18 GeV, \fj\ < 0.7, 
Pt(D*) > 3 GeV and \rj(D*)\ < 1.5, using 
an integrated luminosity of 120 pb~ x . Here Mjj 
is the jet-jet effective mass and fj is the average 
pseudorapidity of the two jets. The fraction of the 
photon momentum producing the two jets is given 

by: x° BS = Sj *2^ e e ~" ■ This quantity is close to 
1 for direct photons and is less than 1 for resolved 
photons. A sample of direct- (resolved-)enriched 
events is defined as x° hs > 0.75 ( < 0.75). 

The differential cross section da/dx^ is com- 
pared in fig. 4 with LO predictions of the MC 
programs PYTHIA Q3|, HERWIG Q3| and CAS- 
CADE and with NLO FO predictions Q]. A 
significant contribution 40%) arises from re- 
solved photons. There is a good agreement in 
shape between the data and the LO MCs, except 
that CASCADE is too high for the high x° bs re- 
gion. The low a;° bs tail of the NLO prediction is 
below the data. 

The angle between the jet-jet axis and the 
beam axis in the dijet rest frame can be approx- 
imated by cos 9* = tanh '' Jetl ~'' 3< " 2 . In QCD the 
differential cross section da/d\ cos#*| is sensitive 
to the spin of the propagator in the hard subpro- 
cess. Whenever the propagator in the LO dia- 
gram is a quark, the cross section rises slowly as 
(1 — | cos 9* |) _1 , and when it is a gluon, it rises 
steeply as (1 - |cos0*|)~ 2 . 

Fig. 5 shows da/d\cos9*\ separately for the 
resolved-enriched (x° hs < 0.75) and direct- 
enriched (x° hs > 0.75) samples. The LO MCs 
PYTHIA and HERWIG describe the shapes of 
both data samples. The direct sample rises slowly 
with | cos 9* | , consistent with the q-exchange 
BGF up-left diagram of fig. 4. The resolved sam- 
ple rises strongly with | cos 9* | , indicating a g- 
exchange signature. Consequently, the LO sub- 
process gg —> cc (up-right diagram of fig. 4) can- 
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not dominate the resolved photon charm-dijet 
process. This suggests the dominance of the LO 
charm-excitation process eg — > eg (down-right di- 
agram of fig. 4). 

The two jets can be distinguished by associat- 
ing the D* meson to the closest jet in rj — <p space. 
Calling this "D* jet" jet 1, the da /d cos 9* rise 
can be studied separately for the photon and pro- 
ton directions. Fig. 7 shows the differential cross 
sections as a function of cosO* for the resolved- 
and direct-enriched samples. The PYTHIA es- 
timation of the resolved process contribution to 
the direct-enriched sample (fig. 7b) explains the 
asymmetric distribution in cos 9* . The strong rise 
in da jd cos 6* towards the photon direction for 
the resolved-enriched sample is clear evidence for 
charm from the photon in the LO picture. The 
NLO FO predictions describe the data well for 
x° hs > 0.75 (fig. 7d), but are below the data 
for x° hs < 0.75 (fig. 7c), both for the proton 
and photon directions. The data shapes are re- 
produced. The CASCADE predictions exceed the 
data by w 30%, mostly for x° bs > 0.75. Again 
the shapes are reasonably well reproduced. 

4. Charm fragmentation 

Charm quark hadronisation into charm mesons 
is parametrised by fragmentation functions, 
which exist in many forms with tunable parame- 
ters fixed from fits to e + e~ data. A direct mea- 
surement of the fragmentation function at HERA 
can reduce theoretical uncertainties and test the 
universality of charm fragmentation. 

Charm-jet PHP events have been measured 
with the ZEUS detector, using an integrated lu- 
minosity of 120 yb~ x , in the kinematic region 
Q 2 < 1 GeV 2 , 130 < W < 280 GeV. At least 
one jet had to satisfy E J T et > 9 GeV, \rf et \ < 2.4 
and a D* meson with pt{D*) > 2 GeV, 
\r](D*)\ < 1.5 had to be associated with a 
jet. The fraction of jet energy carrried by the 
D* was defined as z = where p|| is 

the D* longitudinal momentum relative to the 
jet axis. The normalised differential cross sec- 
tion in z was compared to the PYTHIA LO 
MC with the Peterson fragmentation function ^1 
f(z)cc [z(l — 1/z — e/(l — z)) 2 ] -1 , for various 



values of the free parameter e. Strong sensitiv- 
ity to the e value was found. A fit to the best 
value yielded e = 0.064 ± O.Q06±g;gJ|, compared 
to 0.053 from LO fits to the LEP data. 

The z distribution measured with the ZEUS 
data was compared to that from e + e~ colliders, 
where z was defined as z = Ed* / Ebeam- The 
comparison with OPAL and ARGUS results is 
shown in fig. 6. Similar shapes are obtained for 
z > 0.3 with precision of the HERA data com- 
petitive with the LEP data. The low-z peak in 
the OPAL data is due to gluon splitting to cc. 
The results support the universality of the charm 
fragmentation function. 

The HI measurements jS] of inclusive charm 
mesons in the DIS regime (section 2) were used to 
deduce preliminary fragmentation fractions of the 
charm quark to the various mesons. The results, 
compared with the e + e~ world average values (in 
brackets), are: 

f(c -> D+) = 0.20 =b 0.02±g:g|±g;31 (0.23±0.02) 
f(c -» D°) = 0.66±0.05±g-^t°05 (0.55±0.03) 
f(c - £>+) = 0.16±0.04±°;°ft°;°5 (0.10±0.03) 

f(c -> D*+) = 0.26±0.02±8;g|±g;o2 (0.24±0.01) 
where the first uncertainty is statistical, the sec- 
ond is systematic and the third is from theory. 

In addition, the HI measured cross sections 
were used to test the isospin invariance of the 
fragmentation process by calculating the ratio of 
neutral (cu) to charged (cd) D-meson production 

R _ a dir (D a )+a(D'°) _ a{D Q )-a(D* + )xBR 
n u/d — a d "'(D±)+a(D*±) ~ cr< . D ± )+<j{D' ± ) x BR 

and to extract the strangeness suppression factor 

7s = /( c Jj1 C )+/(c^d") and tnc fraction 
of D mesons produced in a vector state 
P v = vTps = e£td- Here Si? = 0.677 ±0.005 
is the D* + — > D°ir + branching ratio. The results, 
compared with previous values, are: 
R u/d = 1.26 ± 0.20 ± 0.12 (HI prel.) 

1.00 ± 0.09 (e + e~ world average) 

7s = 0.36 ±0.10 ±0.08 (HI prel.) 
0.27 ±0.04 ±0.03 (ZEUS) 
0.26 ± 0.03 (e + e~ world average) 

P v = 0.69 ±0.04 ±0.01 (HI prel.) 

0.55 ±0.04 ±0.03 (ZEUS prel.) 
0.60 ±0.03 (e + e~ world average). 

All the HERA measurements are in good agree- 



ment with e + e world average results, indicating 
again the universality of charm fragmentation. 
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Figure 4. Upper part: LO QCD charm- 
production diagrams. Direct photon, BGF 
(7 g — > cc) ; resolved photon (g g — > cc) ; resolved- 
photon charm excitation (c g — > g c, c in pro- 
ton hemisphere); resolved-photon charm excita- 
tion (c g — > c g, c in photon hemisphere). Lower 
part: Differential cross-section da/dx° hs for the 
data compared with MC simulations (upper plot) 
and NLO FO predictions (lower plot). Each MC 
distribution is normalised to the data, as indi- 
cated in the brackets. The NLO uncertainty after 
hadronisation is given by the shaded band. 
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Figure 5. Differential cross sections da/d\ cos 9*\ 
compared with MC simulations for samples en- 
riched in resolved- (upper plot) and direct- (lower 
plot) photon events. 
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Figure 7. Differential cross sections da /d cos 9* 
compared with MC simulations and NLO FO pre- 
dictions for samples enriched in resolved- (a,c) 
and direct- (b,d) photon events. The shaded ar- 
eas in (a,b) are the contribution of the PYTHIA 
direct photon process to the resolved-enriched 
sample and of the resolved photon process to the 
direct-enriched sample, respectively. The NLO 
uncertainty after hadronisation is given by the 
shaded band. 



Figure 6. Fragmentation function versus z for 
the ZEUS data compared to e + e~ measurements. 
The data sets were normalised to 1/ (bin width) 
for z > 0.3. 



